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SUMMARY 


This  memorandum  addresses  the  problem  of  locating  linear  features,  not 
necessarily  straight,  in  Synthetic  Aperture  Radar  images.  The  characteristics 
of  such  features  and  the  way  in  which  these  characteristics  affect  the  desired 
properties  of  line  finding  algorithms  are  discussed. 

Two  line  detection  algorithms  are  then  described:  the  Hough  transform, 
and  the  dynamic  programming  technique.  Results  of  applying  both  algorithms 
to  simulated  and  real  SAR  images  are  given.  The  Hough  transform  is  found  to 
be  too  restricted  to  be  useful  for  this  application.  The  dynamic  programming 
technique  locates  lines  well.  The  problem  of  spurious  lines  appearing  in 
clutter  is  discussed. 
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INTRODUCTION 


Remotely  sensed  images  of  terrain  contain  a  large  number  of  different 
features:  fields,  trees,  hedges,  rivers,  buildings,  roads  etc.  The 

identification,  extraction  and  analysis  of  interesting  features  from  images 
can  be  a  very  time  consuming  task  if  done  manually.  If  the  image  sensor 
produces  data  at  a  significant  rate,  this  task  may  be  prohibitively  expensive. 

The  RSRE  X-band  SAR  is  one  example  of  such  a  system:  images  can  be  produced 
in  real  time  over  a  swath  10  km  wide,  at  an  aircraft  velocity  of  the  order  of 
200  m/s.  The  resolution  of  the  system  is  of  the  order  of  metres,  and  thus 
images  are  produced  at  the  rate  of  approximately  one  million  pixels  per 
second.  Interpretation  of  the  images  is  complicated  still  further  by  the  presence 
of  speckle  noise  which  is  a  consequence  of  the  coherent  nature  of  the  imagine 
system. 

Since  human  interpretation  of  the  complete  set  of  image  data  is  not 
feasible,  there  is  a  requirement  for  an  automatic  system  of  image  analysis. 

Even  if  this  system  is  unable  to  perform  a  complete  analysis  of  the  image,  it 
may  still  be  used  to  select  a  subset  of  the  image  that  can  be  presented  to  a 
higher  level  operator  (either  a  human,  or  a  more  intelligent  system  that  is 
too  slow  to  be  used  for  analysis  of  the  complete  image). 

This  memorandum  describes  some  initial  work  that  has  been  done  towards 
defining  one  part  of  the  low  level  image  analysis  system:  namely  that  of 
detecting  linear  features  in  images.  Examples  of  such  features  are  roads, 
hedges  and  fences.  Road  detection  is  a  worthwhile  objective,  since  it  is 
likely  that  other  objects  of  interest  may  be  found  in  their  vicinity,  particularly 
vehicles.  Also,  linear  features  can  be  used  for  image  to  image  and  image  to  map 
registration. 

Section  2  discusses  the  effect  of  the  speckle  noise  in  SAR  images  on  the 
appearance  of  lines  in  images,  and  the  way  this  influences  the  required 
properties  of  a  line  finding  algorithm.  Sections  3  and  4  discuss  two  particular 
line  detection  algorithms  in  detail:  the  Hough  transform  and  dynamic  programming 
respectively.  The  results  of  applying  the  two  algorithms  to  synthetic  and  real 
image  data  are  presented  in  Section  5.  Conclusions  and  suggestions  for  further 
work  are  contained  in  Section  6. 

2  THE  EFFECTS  OF  SPECKLE 


SAR  is  a  coherent  imaging  system.  The  major  consequence  of  this  fact  is 
that  the  images  are  degraded  by  speckle  noise.  This  noise  is  uncorrelated  from 
resolution  cell  to  resolution  cell,  and  the  probability  density  function  of  the 
intensity  over  a  region  of  constant  radar  cross  section  is  given  by 


P(I)dI  =  <T>  exp(- dl 


(2.1) 


where  <I>  is  the  mean  intensity  over  the  region. 


Suppose  we  have  a  region  with  mean  intensity  <It>,  against  a  background 
of  mean  intensity  <Ib>.  Then  we  can  set  a  threshold  for  detection  of  the  region 
against  the  background  at  the  point  at  which  the  two  probability  distributions 
are  equal,  given  by 
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The  probability  that  a  single  pixel  in  the  region  will  be  incorrectly 
classified  as  being  part  of  the  background  is  then  given  by 


D  =  P(IT  <  T) 


(2.3) 


For  example,  for  a  region  whose  mean  intensity  is  10  times  that  of  the 
background,  the  probability  that  a  single  pixel  in  the  region  has  an  intensity 
less  than  the  threshold  is  approximately  0.23.  Fig  2.1  shows  the  behaviour  of 
D  as  a  function  of  (it)/(IbV 


Speckle  noise  also  affects  the  appearance  of  linear  features  in  the 
image.  Consider  a  line  of  mean  intensity  <Il>  against  a  background  of  mean 
intensity  <Ig>.  Then  we  can  define  a  threshold  for  detection  of  the  line  as 
above,  and  the  probability  of  any  single  pixel  on  the  line  having  an  intensity 
below  this  threshold  is  given  by  2.3  above.  We  also  wish  to  calculate  the 
probability  Pr,  of  the  gap  of  length  R  pixels  in  a  line.  Clearly  this 
probability  increases  with  the  length,  N,  of  the  line,  but  it  is  not  possible 
to  find  an  exact  formula  for  Pr. 


In  the  absence  of  an  analytic  formula  for  Pr,  a  program  was  written  to 
investigate  the  behaviour  of  Pr,  as  D,  R  and  N  were  varied.  This  was  done  by 
generating  sequences  of  up  to  N  binary  numbers  with  probability  D  of  being  set. 
If  a  sequence  of  R  set  numbers  was  detected,  a  success  was  recorded.  10000 
trials  were  performed,  and  the  probability  of  a  gap  of  length  R  in  a  line  of 
length  N  was  thus  given  by  the  number  of  successes  divided  by  10000.  5  values 
of  D  were  used,  and  the  probabilities  of  gaps  of  length  1,  2,  3,  h  and  5  pixels 

in  lines  of  lengths  5,  10,  15 . 50  pixels  were  estimated.  These  results  are 

tabulated  in  appendix  1,  and  Fig  2.2  shows  the  probabilities  of  gaps  of 
various  lengths  for  different  length  lines  when  the  probability  of  a  single 
pixel  "dropping  out"  is  0.2,  corresponding  to  a  line  which  is  ~12.22  times 
brighter  than  the  background.  Clearly  there  is  a  significant  probability  of  a 
gap  of  length  3  pixels  for  lines  which  are  20  pixels  or  more  in  length. 

These  results  shown  that  a  local  operator  with  a  window  size  of  less  than 
seven  pixels  is  unlikely  to  be  successful  in  locating  lines  on  SAR  images.  An 
algorithm  is  required  that  is  capable  of  taking  a  more  "global"  view  of  the 
image,  and  spanning  the  gaps  in  the  line  due  to  speckle.  Two  such  algorithms 
will  be  discussed  in  the  following  two  sections. 


3  THE  HOUGH  TRANSFORM 

The  Hough  transform  (eg  Duda  and  Hart  1972)  is  a  widely  used  line 
detection  algorithm.  It  is  a  global  procedure,  in  that  all  points  that  lie  on 
a  given  line  in  an  image  map  onto  a  single  point  in  a  transformed  space. 

3.1  IMPLEMENTATION 

Any  straight  line  can  be  parameter ised  by  the  angle  that  the  normal 
to  it  makes  with  the  horizontal,  0,  and  its  distance  r,  from  an  origin, 
typically  the  centre  of  the  image.  (The  "standard"  parameterisat ion  of 
a  line,  in  terms  of  its  slope  and  intercept  with  an  axis  is  not 
appropriate  here,  since  these  quantities  can  increase  without  limit). 

The  Hough  transform  proceeds  as  follows:  each  point  in  the  image 
lies  on  many  possible  straight  lines,  and  thus  maps  onto  many  points  in 
r,  0  space.  For  an  image  point  at  x,y,  the  corresponding  points  in  the 
transform  space  are  given  by 
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r  =  x  cos(0)  +  y  sin  (0) 


(3.1) 


i.e,a  sinusoidal  curve.  The  image  intensity  at  (x,y)  is  added  to  the 
transform  for  each  point  on  the  curve.  This  procedure  is  repeated  for 
all  points  in  the  image.  If  many  points  lie  on  the  same  straight  line, 
their  corresponding  curves  will  intersect  at  a  particular  point,  leading 
to  a  peak  in  the  Hough  transform.  The  location  of  these  peaks  gives 
the  parameters  (r,0)  of  the  lines  in  the  image. 

The  fact  that  all  points  in  the  image  are  used  independently  makes 
the  method  robust  to  missing  points  in  the  line  -  their  effect  is  simply 
to  reduce  the  height  of  the  peak  in  the  transform  space.  The  discussion 
of  section  2  showed  that  this  is  a  desirable  quality  for  a  SAR  image 
line  finding  algorithm. 

When  locating  peaks,  it  is  important  to  take  into  account  the  fact 
that,  depending  on  the  angle  and  distance  of  the  line  from  the  origin, 
different  lengths  of  the  line  intersect  the  image,  and  so  make  different 
contributions  to  the  Hough  transform.  For  example,  a  line  passing 
diagonally  through  the  centre  of  the  image  clearly  has  a  much  longer 
intersection  than  one  which  passes  through  one  corner.  For  this  reason, 
the  Hough  transform  must  be  normalised  by  dividing  it  by  the  transform 
of  a  uniform  image  before  being  analysed. 

The  Hough  transform  space  is  quantised:  if  more  samples  are  taken 
in  0,  then  the  orientation  of  lines  can  be  found  more  accurately,  while 
more  samples  in  r  gives  more  accurate  positions.  However,  the  more 
samples  there  are  in  the  Hough  transform,  the  longer  it  will  take  to 
compute:  for  an  image  that  is  NX  by  NY  pixels,  with  N0  samples  in  0, 

and  NR  samples  for  r,  the  number  of  operations  required  is  proportional 
to  NX. NY. N0. NR.  For  initial  processing  of  an  image,  coarse  sampling  in 
r  and  0  is  acceptable  in  order  to  find  approximate  lines.  Finer  sampling 
can  then  be  used  if  more  accurate  values  for  the  parameters  of  particular 
lines  are  required. 

Figure  3.1  shows  an  image  consisting  of  a  single  straight  line  and 
Fig  3.2  its  normalised  Hough  transform.  The  peak  in  the  transform 
corresponding  to  a  line  having  r=20,  and  0=tt/6  is  clearly  visible. 

Having  obtained  the  Hough  transform,  we  then  have  to  "invert"  it  to 
obtain  the  lines  in  the  image.  However,  it  is  not  possible  to  do  this 
directly,  since  the  location  of  a  peak  in  the  Hough  transform  gives  the 
orientation  and  position  of  the  line,  but  contains  no  information  about  its 
end  points.  To  overcome  this  problem,  the  following  algorithm  was 
developed.  For  a  given  line,  all  the  pixels  lying  on  this  line  in  the  image 
were  extracted  into  a  single  vector.  Then,  for  all  possible  start  and 
end  points,  subject  to  the  line  being  a  minimum  length,  in  this  vector,  the 
mean  intensity  was  calculated.  The  positions  of  the  start  and  end  points 
corresponding  to  the  maximum  mean  were  recorded,  and  used  as  the  start  and 
end  points  for  this  line  segment.  When  a  pixel  was  assigned  to  a  line, 
its  contribution  to  the  Hough  transform  was  removed,  so  that  further 
peaks  in  the  transform,  corresponding  to  different  lines,  could  be 
located . 
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DYNAMIC  PROGRAMMING 


Section  2  showed  that  a  successful  line  finding  algorithm  must  be 
capable  of  spanning  gaps  of  approximately  three  pixels.  To  do  this,  it  must 
take  a  "global"  view  of  the  image.  One  way  in  which  this  could  be  done  in 
principle  is  to  generate  all  possible  line  segments  of  length  N  pixels  in  the 
image,  and  to  calculate  a  figure  of  merit  for  each  segment.  The  figure  of 
merit  could  be  the  summed  intensity  along  the  line,  with  a  factor  that 
penalises  highly  curved  1 ines .  However,  this  approach  is  not  feasible,  since 
even  for  an  image  of  only  64  by  64  pixels,  there  are  clearly  far  too  many 
possible  line  segments  of  length  10  pixels,  say,  for  the  figure  of  merit  to 
be  calculated  for  all  of  them.  The  method  of  dynamic  programming  (eg  Montanari 
1971),  overcomes  this  combinatorial  explosion. 

4.1  IMPLEMENTATION 

Suppose  the  figure  of  merit  for  a  line  seg  lent  of  length  N  is  the 
summed  intensity  over  the  line,  with  a  penalty  for  highly  curved  lines. 

Then  the  desired  line  segment  is  given  by  the  values  of  zj,  Z2, . 

and  d^,  d2> . d^  which  maximise 

N_^ 

I(z.)  -  q  /  (d.-d.  Jmod  8  (4.1) 

J  J  J~> 

j-2 

where  zj=(xj,yj)  is  the  coordinate  in  the  image 

I(zj)  is  the  intensity  at  pixel  (xj,  y j ) , 
q  is  the  penalty  factor  for  curves  in  the  line, 
dj  is  the  exit  direction  from  the  pixel  zj,  which  lies  in  the 
range  0  to  7,  defining  the  directions  as  follows 

0 

7  1 

6  z;  2 
5  3 

4 

A  pair  of  variables  zj,  dj  are  referred  to  as  a  state. 

There  are  two  constraints  that  must  be  satisfied  by  the  variables  zj 
and  dj: 

max(|xj  -  x j  —  1 1  , | y j  -  yj-l|)  =  li  which  ensures  that  pixels  on  a  line 

are  contiguous. 

((dj  -  d j_ j )mod8 )  < 1 ,  which  prevents  the  line  from  turning  through 

a  right  angle  in  one  pixel. 

The  calculation  proceeds  over  N  stages  as  follows.  At  stage  j,  we 
have  to  decide  whether  a  state  zj,  dj  should  be  associated  with  the  jth 
point  of  the  curve.  We  make  this  decision  on  the  basis  of  the  profit  of 
associating  this  state  with  the  (j-l)**1  point  of  the  curve.  Thus  at  each 
stage  we  compute 

Vj  (z^,  d|^)=max(  I(z^)  -  q  (d^  -  dj)mod8)  +  V  j_  ^  ( z^ ,  d  j ) , 


where 


dfc  is  the  exit  direction  from  z{, 

d]  is  the  entry  direction  into  z£,  and  the  exit  direction 
from  Zh, 

Zh  is  the  neighbour  of  along  dj. 

The  calculation  starts  with  V0(z£,.)  =  Kz^).  At  each  stage,  a 
memory  is  updated  which  records  the  path  to  a  given  state.  This  is  used 
to  trace  back  along  the  line. 

Having  completed  N  stages,  we  find  the  largest  profit  value,  which 
gives  the  end  point  of  the  best  line  in  the  image: 

PN  =  xi>yi  st-  V*i'  ’)=max 

zi 

Having  found  the  end  point  of  this  line,  we  can  then  trace  back  through 
the  states  that  brought  us  to  this  point,  using  the  memory  described  above. 

At  each  stage  in  the  traceback,  when  a  pixel  is  included  in  a  line,  its 
effect  on  the  profit  memory,  V,  is  removed,  so  that  it  cannot  be  included 
in  another  line. 

This  algorithm  is  expressed  in  a  stylised  language  in  appendix  2. 

It  will  be  clear  that  the  method  makes  rather  large  demands  on  memory  to 
record  the  traceback  information,  and  also  requires  a  large  number  of 
operations.  To  find  lines  of  length  N  in  an  image  of  MX  by  MY  pixels 
requires  N. MX. MY.  8. 3  traverses  of  the  inner  loop.  The  factors  8  and  3 
arise  from  the  fact  that  there  are  8  possible  exit  directions  from  a 
pixel,  and  for  each  exit  direction  there  are  3  possible  entry  directions. 

The  inner  loop  requires  one  addition  and  one  subtraction  of  image  values, 
which  can  be  stored  as  fixed  point  numbers.  Bertolazzi  and  Pirozzi  (1984) 
describe  how  the  algorithm  can  be  implemented  on  a  parallel  architecture: 
if  there  are  sufficient  processors  such  that  one  pixel  can  be  associated 
with  one  processor,  then  the  algorithm  runs  in  time  0(N).  If  there  are 
P  processors,  the  run  time  is  0(N.MX.MY/P) ,  ie  the  algorithm  increases 
in  speed  proportional  to  the  number  of  processors  available.  The  memory 
requirement  of  the  algorithm  is  mainly  determined  by  the  traceback  memory. 
Although  the  profit  values  must  be  recorded,  they  do  not  need  to  be  stored 
for  every  stage,  but  can  be  double-buffered  so  that  the  profit  at  the 
current  stage  can  be  evaluated  from  the  previous  stage.  This  double 
buffer  uses  2.  MX. MY. 8  storage  locations.  For  image  values  in  the  range 
0-255  (ie  8  bit  data),  and  for  N=16,  for  example,  each  storage  location 
will  require  12  bits.  The  memory  required  for  traceback  needs  to  be 
MX.  MY.  8. N  locations,  but  each  location  need  be  only  three  bits,  since  it 
holds  a  direction  in  the  range  0-7. 

5  RESULTS 

This  section  presents  the  results  obtained  when  the  algorithms  described  in 
sections  3  and  4  were  applied  to  synthetic  and  real  image  data.  For  both  types 
of  data,  a  threshold  was  applied  to  the  images  as  follows:  the  mean  value  of 
the  image  was  calculated,  and  only  points  that  exceeded  this  mean  value  were 
used  in  calculations.  This  has  the  effect  of  reducing  the  run  time  of  the 
programs . 

5.1  SYNTHETIC  DATA 

The  test  image  used  in  this  case  consisted  of  a  spiral  imposed  on  a 
background.  The  intensity  of  the  spiral  was  10  times  that  of  the  background. 
This  form  of  test  image  was  chosen  as  it  contains  line  segments  of  differing 
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curvatures.  Fig  5.1  shows  the  test  image,  displayed  in  modulus  form  to 
reduce  the  dynamic  range  of  the  image.  A  thresholded  version  of  the 
image  is  shown  in  Fig  5.2  with  the  threshold  set  according  to  the  results 
in  section  2.  This  shows  the  dropouts  that  occur  in  the  line. 

5.1.1  HOUGH  TRANSFORM  APPLIED  TO  SYNTHETIC  DATA 

The  Hough  transform  of  the  image  is  shown  in  Fig  5.3.  The 
stripes  in  the  transform  correspond  to  lines  in  the  image  whose 
distance  from  the  origin  increases  as  the  angle  increases,  as  would 
be  expected  for  a  spiral.  The  reconstructed  image  is  shown  in  Fig  5.4. 
All  points  in  the  Hough  transform  having  a  value  greater  than  the  mean 
of  the  original  image  were  used  to  generate  the  lines  shown  in  the 
reconstructed  image.  This  corresponds  to  reconstructing  lines  whose 
mean  intensity  is  greater  than  the  image  mean  value. 

The  line  segments  towards  the  edge  of  the  spiral  have  been 
located  fairly  successfully,  while  the  more  central  ones  have  not. 

This  is  to  be  expected,  since  the  Hough  transform  is  matched  to  straight 
lines,  and  it  is  only  at  the  edge  of  the  spiral  that  the  line  segments 
appear  reasonably  straight.  This  is  a  major  problem  with  the  Hough 
transform  as  implemented  here;  it  can  only  locate  straight  line  segments. 
It  is  possible  to  modify  the  Hough  transform  to  locate  any  parameteris- 
able  shape,  but  it  cannot  find  a  line  of  arbitrary,  and  possibly 
variable  curvature. 

5.1.2  DYNAMIC  PROGRAMMING  APPLIED  TO  SYNTHETIC  DATA 

We  now  examine  the  results  of  applying  the  dynamic  programming 
algorithm  to  the  synthetic  image.  As  described  in  section  4,  this 
method  is  capable  of  locating  curved  lines.  Fig  5.5  shows  the  lines 
located  in  the  synthetic  image  when  no  curvature  penalty  was  applied 
(ie  q=0).  Although  some  line  segments  have  been  located  correctly, 
there  are  many  loops  where  the  algorithm  has  joined  a  pixel  to  itself 
via  the  shortest  possible  route.  To  eliminate  these  loops,  the  curvature 
penalty  was  applied,  with  q=50  (the  image  values  are  scaled  in  the 
range  0-255).  The  reconstructed  image  obtained  in  this  case  is  shown 
in  Fig  5.6.  The  loops  have  beem  eliminated  using  the  curvature  penalty, 
although  it  must  be  remembered  that  the  curvature  penalty  wil1  prevent 
the  algorithm  locating  highly  curved  lines. 

5.2  SAR  IMAGE  DATA 

5.2.1  PREPROCESSING  OF  RSRE  SAR  IMAGES 

Fig  5.7  shows  a  section  from  a  SAR  image  of  a  road  at  a  nominal 
resolution  of  a  few  metres  in  range  and  azimuth.  The  hedges  on  each 
side  of  the  road  are  clearly  visible.  It  will  be  noticed  that  there 
are  correlations  in  the  vertical  direction  in  the  image,  corresponding 
to  lines  of  constant  azimuth  position.  These  are  due  to  the  fact 
that  the  range  dimension  sampling  is  about  two  thirds  of  the  range 
resolution.  In  order  to  remove  these  range  correlations,  the  frequency 
spectrum  of  the  image  in  the  range  direction  was  whitened,  as  follows: 
a  large  image  was  selected  which  contained  relatively  little  structure. 
From  this  image,  vertical  lines  were  extracted.  Each  line  was 
Fourier  transformed,  and  the  resulting  spectra  were  summed  in  power. 

The  spectrum  of  the  image  in  Fig  5.7  was  then  obtained,  and  each  line 
of  constant  azimuth  frequency  was  divided  by  the  root  mean  square 
power  spectrum  obtained  from  the  large  image.  This  modified  Fourier 
transform  was  then  inverted  to  give  the  decorrelated  image.  The  result 


of  this  process  is  shown  in  Fig  5.8.  It  is  clear  that  the  range 
correlations  have  been  significantly  reduced  while  preserving  the 
structure  of  the  image.  The  line  finding  algorithms  were  then 
applied  to  this  image. 

5.2.2  HOUGH  TRANSFORM  APPLIED  TO  SAR  IMAGE  DATA 

The  Hough  transform  of  Fig  5.8  is  shown  in  Fig  5.9.  The  peaks 
in  the  transform  due  to  the  hedges  on  each  side  of  the  road  are 
clearly  visible.  The  reconstructed  image  from  this  transform  is 
shown  in  Fig  5.10.  Although  portions  of  the  hedges  have  been  located, 
the  majority  have  been  missed.  This  is  because  the  road  is  curved, 
resulting  in  broadened,  lower  peaks  in  the  Hough  transform  space. 

This  again  demonstrates  the  "tuned"  nature  of  the  Hough  transform. 


5.2.3  DYNAMIC  PROGRAMMING  APPLIED  TO  SAR  IMAGE  DATA 
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The  problem  of  lines  being  located  in  clutter  also  raises  the 
question  of  when  the  traceback  process  should  be  terminated.  The 
traceback  is  started  from  the  pixel  having  the  highest  current  value 
in  the  profit  array,  V.  Apart  from  the  curve  penalty  factor,  q,  this 
value  divided  by  the  number  of  stages  is  the  mean  value  along  the 
current  line.  This  gives  an  indication  of  the  confidence  in  the  line, 
and  its  strength  relative  to  the  background,  and  hence  can  be  used  to 
determine  whether  the  traceback  should  be  done  for  the  line.  However, 
a  single  very  bright  point  in  clutter,  due  to  an  isolated  strong 
target  for  example,  can  give  rise  to  reasonable  mean  values  over  a 
line.  For  this  reason,  a  further  "line  editing"  step  is  probably 
required,  in  which  each  line  segment  located  by  the  algorithm  is 
analysed  to  ensure  that  it  has  reasonable  properties.  One  simple 
test  is  to  measure  the  variance  along  the  line;  if  this  is  large,  then 
the  line  contains  a  few  bright  pixels,  and  should  probably  be  rejected. 
More  complicated  tests  could  also  be  devised  based  on  the  context  of 
the  line  image.  For  example,  if  the  line  is  thought  to  be  a  hedge, 
then  a  shadow  would  be  expected  behind  it. 

6  CONCLUSIONS 

Two  line  location  algorithms  have  been  studied  in  this  memorandum.  The 
Hough  transform  was  found  to  locate  straight  lines  well.  However,  the  fact  that 
interesting  lines  in  SAR  images  need  not  be  straight  limits  the  usefulness  of 
the  Hough  transform.  The  dynamic  programming  technique  locates  curved  and 
straight  line  segments  well,  but  is  expensive  in  terms  of  storage  and  epu  time 
required.  A  multiprocessor  implementation  of  the  technique  can  reduce  the  run 
time  significantly. 
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Apart  from  the  problems  of  gaps  in  lines,  speckle  can  also  give  rise  to 
spurious  lines  in  the  background.  These  lines  can  be  suppressed  by  demanding 
that  the  lines  located  have  a  minimum  length.  For  the  lines  segments  studied 
here,  which  had  a  target  to  background  intensity  ratio  of  approximately  10  to  1, 
it  was  found  that  a  line  had  to  be  approximately  20  pixels  long  before  it  could 
be  located.  Further  work  could  investigate  an  "editing"  process  to  be  applied 
to  the  lines  to  eliminate  these  spurious  lines. 
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APPENDIX  1.  PROBABILITY  OF  GAPS  FOR  LINES  OF  LENGTH  N 
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APPENDIX  2 


This  appendix  presents  the  implementation  of  the  dynamic  programming 
algorithm  in  a  stylised  language. 

for  Each  Stage 

for  each  pixel 

for  each  exit  direction  from  this  pixel 

for  each  entry  direction  for  this  exit  direction 
if  [(profit  from  previous  stage 
+  image  value  for  this  pixel 

-  curvature  penalty  for  given  entry  and  exit  direction) 
is  a  maximum]  then 

record  the  maximum  profit 

record  the  entry  direction  for  this  exit  direction 

end  if 

next  entry  direction 
next  exit  direction 
next  pixel 

next  Stage 

set  threshold  for  reconstruction  of  lines 
loop: 

scan  the  profit  array  for  its  maximum  value 
if  (maximum  profit  exceeds  threshold]  then 

record  the  pixel  corresponding  to  the  maximum  profit 

obtain  the  direction  line  enters  this  pixel  from  the  memory 

remove  the  effect  of  this  pixel  from  the  profit  array 

for  each  stage  starting  from  the  last 


get  the  next  pixel  on  the  line  from  the  previous  pixel  and  its  direction 
obtain  the  direction  line  enters  this  pixel  from  memory 
trace  the  effect  of  this  pixel  back  to  the  final  stage 


remove  its  effect  from  the  profit  array 
next  stage 

else 

8  top 

end  if 
goto  loop 
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reconstruction  of  Fig.  5.8 
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Abstract 

This  memorandum  addresses  the  problem  of  locating  linear  features,  not 
necessarily  straight,  in  Synthetic  Aperture  Radar  images.  The  characteristics  of 
such  features  and  the  way  in  which  these  characteristics  affect  the  desired 
properties  of  line  finding  algorithms  are  discussed. 

Two  line  detection  algorithms  are  then  described:  the  Hough  transform,  and  the 
dynamic  programming  technique.  Results  of  applying  both  algorithms  to  simulated 
and  real  SAR  images  are  given.  The  Hough  transform  is  found  to  be  too  restricted 
to  be  useful  for  this  application.  The  dynamic  programming  technique  locates 
lines  well.  The  problem  of  spurious  lines  appearing  in  clutter  is  discussed. 


